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Abstract 

We consider the Euler equations of incompressible inviscid fluid dynamics. We 
discuss a variational formulation of the governing equations in Lagrangian co¬ 
ordinates. We compute variational symmetries of the action functional and 
generate corresponding conservation laws in Lagrangian coordinates. We clar¬ 
ify and demonstrate relationships between symmetries and the classical balance 
laws of energy, linear momentum, center of mass, angular momentum, and the 
statement of vorticity advection. Using a newly obtained scaling symmetry, we 
obtain a new conservation law for the Euler equations in Lagrangian coordinates 
in n-dimensional space. The resulting integral balance relates the total kinetic 
energy to a new integral quantity defined in Lagrangian coordinates. This re¬ 
lationship implies an inequality which describes the radial deformation of the 
fluid, and shows the non-existence of time-periodic solutions with nonzero, finite 
energy. 
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1. Introduction 


This study concerns the Euler equations of incompressible inviscid fluid dy¬ 
namics. The Euler equations provide an accurate representation of a variety of 
5 inviscid fluid flows and are used in numerous practical situations. Mathematical 
interest in the Euler equations includes the open problems of classifying nonlin¬ 
ear blow-up, long-time behavior, and stability of solutions for various types of 
initial conditions. Conservation of energy plays an important role in controlling 
blow-up; for example, solutions of infinite energy have been shown to blow-up 
10 in a finite amount of time [l| . 

The goal of this study is to use variational symmetries of the Euler equations 
in Lagrangian coordinates to find conservation laws. Noether’s First Theorem 
Q relates variational symmetries to conservation laws of a corresponding dif¬ 
ferential system. However, the Euler equations in Eulerian variables do not 
15 admit a variational formulation. Point symmetries of the Euler equations in 
Eulerian variables were obtained in [^; see also P|. There are many works 
on the conservation laws of the Euler equations; see Q] and references therein; 
see also [g, 7, 8| for conservation laws in n-space computed using the direct 
method. Conservation laws were obtained for the constrained Euler equations 
20 in Q, and for the Euler equations in vorticity formulation 


hand, the Euler equations in Lagrangian variables (see 
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do admit a variational formulation. Caviglia and Morro 


. 101. On the other 

IT 

12| for background) 


13j | computed symme¬ 


tries and conservation laws using a variational formulation for the compressible 
Euler equations in Lagrangian coordinates. However, they later [l^ computed 
conservation laws for the incompressible Euler equations without using sym¬ 
metries or a variational formulation. We show that a simple modification to 
their compressible flow Lagrangian in [l^ allows for a variational formulation 
of incompressible flow. 

In this paper we study a variational formulation of the incompressible Euler 
equations in Lagrangian coordinates. We obtain a new scaling symmetry of 
the action functional. This scaling symmetry leads via Noether’s theorem to a 
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new conservation law in n-space. The resulting conserved quantity relates the 
total energy integral, a quantity defined in Eulerian variables, to a new integral 
quantity defined exclusively in Lagrangian variables. This relationship allows 
35 for the quantification of the radial deformation of the fluid, and rules out the 
existence of time-periodic solutions with nonzero, finite energy. 


2. Variational Framework 

We first introduce Lagrangian coordinates, discuss their physical interpre¬ 
tation, and detail explicitly the transformation of the Euler equations to La- 
40 grangian coordinates from Eulerian variables. Then, we demonstrate the varia¬ 
tional formulation behind the Euler equations in Lagrangian coordinates using 
a Lagrange multiplier (namely the pressure) of the incompressibility constraint. 


2.1. Lagrangian Coordinates 


The contin uity and Euler equations of incompressible ideal fluid dynamics 
are as follows 15j: 


Va; • u = 0, 

[dt +u-Va:)u = -VxP, 


( 1 ) 

( 2 ) 


where u = (m^, ..., u") is the velocity vector, Vx = {dx ^, • ■ •, dx^) is the gra- 
45 dient, and p = p/p is the pressure p divided by the constant density p; we 
henceforth refer to p as the pressure. The velocity vector u and pressure p are 
functions of {x,t) = (a;^,..., x”, t). The first condition ([T]) gives the local con¬ 
servation of mass and specifies that the fluid be incompressible. The second 
states the momentum balance within a small parcel of fluid. The equations 
50 are considered on n -|- 1 dimensional spacetime, where n > 2 is 2 or 3 in most 
applications. 

The equations for ideal incompressible fluid dynamics can be recast using the 
“Lagrangian map” x(a, t) = (a;^,..., cc") 1^ l3| ■ This smooth function gives the 
position a: of a particular fluid particle given by “label vector” a = (a^,..., a”) 
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at time t. At the initial time t = 0, the positions of the fluid particles coincide 
with their “Lagrangian coordinates”: a;(d,0) = a. As time progresses, the 
positions of the particles deviate from their initial positions by the action of 
the fluid velocity u{x,t). This relationship between position and velocity can 
be quantified by a kinematic system of equations: dtx{d,t) = u{x,t), where dt 
denotes partial differentiation with respect to variable t. For small times t —>■ 0, 
Taylor’s theorem gives 


x{d, t) « x{d, 0) + t dtx(d, 0) + 0(t^), 

= a + tu{x{d,0),0) + 0{t^) (3) 

= d + t u(d, 0) + 


In this moving frame that follows the fluid particles, governing equations 
ID), ID take the following form [ll]. [l2|: 


d'^x^ j dp da^ 


( 4 ) 

i = (5) 


It is convenient to introduce some notation. Subscripts denote partial deriva¬ 
tives with respect to the fluid labels a® (f = 1 ,..., n) or time t: fj = ^ and 
ft = ^ (superscript indices denote vector or matrix components). Here, = 
55 ^ is the inverse of the Jacobian matrix Xj such that x^^x\. = 

(identity matrix). We assume summation over repeated indices. 

Equation ID can be verified to be equivalent to |D using Jacobi’s formula: 


ddet A = det A tr(A ^dA), 


( 6 ) 


where dA is the infinitesimal change in matrix A, and tr(A*'^ ) = A“ is the trace 
of A. In particular, time differentiating |D gives 

dtdet{x]) = det(x})tr = x^^x^ = (7) 

This expression is zero upon differentiating |D, which recovers |D- That the 
Jacobian is unity follows from its constancy and the initial condition x{a, 0) = a. 
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The left side of equation ([2]) can be recovered from that of ([5]) by using the 
definition of the velocity field and the chain rule for total time derivatives: 


xlt = = 


= ul + 


( 8 ) 


2.2. Extremals of the Action Functional 

We show that Eqs. (|3]) and (O are the variational Euler-Lagrange equations 
for an “action” functional <S given in the following form: 


5[a:,p] = J L{dtx)+pG^x^ dtd"a = J £(^dtx,Vx,pJ 


dfd" 


(9) 


Here, L(dtx) = ^\dtx\^ = \xlx\ is the kinetic energy density, and p{a,t) func¬ 
tions as a Lagrange multiplier (see [f^ l for non-holonomic constraint function 
G = det (cc*) — 1 that accounts for the flow incompressibility (in this con¬ 

text, V denotes a “Lagrangian gradient” V = (S^i,..., 9o"))- The integration 
is taken over a region i? of n -I- 1 dimensional spacetime. In what follows, we let 
Di denote total differentiation with respect to the ith variable. 

A variational formulation has been known for incompressible flow since La¬ 


grange 


H, 


but the action in that case does not involve the pressure as an 


explicit dynamical variable. Instead, the action is taken to be 


5[x] = f L{dtx)dt(r‘a (10) 

JR 

subject to a local incompressibility constraint G{Vx) = 0 for all {x,t) € R; the 
65 pressure p arises indirectly once an extremum of the action is sought. While 
these two formulations yield equivalent solutions to the variational problem, we 
find that performing computations with (jO]) is more transparent. 

Lagrangian £ in ([9]) is a slight modification of the function used in 
C{x,p) = pqL — pqE(p) -|-p(p)(det(a:*) — po/p), where E{p) is the internal 
70 energy, p(a,t) is the density, po(d) = p(d) 0) is the initial density, and p = p{p) 
is the pressure. This is a Lagrangian for the compressible and isentropic Euler 
equations in Lagrangian coordinates. For these equations, the density p is a 
dynamical variable and the pressure p = p{p) is a known function of p. This 
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is in contrast to the incompressible equations we study, for which p = po is a 
75 constant, and p(d, t) is a dynamical variable. One can obtain our Lagrangian in 
(P|) from that in 13| by dividing by the constant p = po, discarding the internal 
energy constant E, and replacing p with p as the dynamical variable that C 
depends on. Note that the same authors in the later paper [l^ claimed there 
was no variational formulation for incompressible flow. 

80 For Eulerian system 0-0, there is no variational formulation. However, 


H 


one can be found for the related system of compressible flow; see ua for this 
formulation as well as another for irrotational compressible flow, where velocity 
u = V(j) is restricted to a potential representation. 

We seek extremals x and p of the action functional S. To this end, we take 
X ^ X + e ff(d, t) and p ^ p + eip{d,t) for e —>■ 0. Arbitrary functions fj and ^ 
are assumed to vanish on dR. The first variation SS of action S is: 
diS = De <S[x + ep, p + etp] I 


le=0 


IR 


dtd" 


illLxi + i^G+pdet (^Vfj x^^pi 
We computed the last term using 

D^G {x\ + ep]) = D^det [x] + epj) = det (^Vx) tr 


( 11 ) 


= det 


(w) 




Since is an arbitrary function, it follows from (Hip that 

G{Xj) = det(x®) —1 = 0, which recovers incompressibility condition (jT]). Letting 
px^^pf = D^{pi^^pi) -pp^xi -p^x^^p^ and plL^i^ = Dt{p'L^i^ - p^DtL^i, Eq. 
CH) becomes: 


6S = [ [Dtip^xD + Di{px^^p^)] dtd^a 


( 12 ) 


The first integral contributes via Gauss’ Theorem to surface terms that vanish 
on OR. Hence, for the first variation SS in (HH) to vanish for arbitrary f]{d,t), 
it follows that 

xl^ + PiX^^ + pxf = 0 
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. To recover ([S]), we show that 

if=0, j = (13) 

when Q holds. Since xj. = differentiation and summation with respect 
85 to i gives xl^ x^ + xj/^ = 0, from which we see that = —i®-^ xj^x'^^, which 
re-expresses the last term in But the gth partial derivative of det(a;®) = 1 
is dq det(xj) = 0 = tr ^x^^xl^"j = i®-^ xj^, so xf^ = —{x^^xig)x'^^ = 0. 

3. Symmetries of the Action 

In this section, we find all (variational) point symmetries of the action func- 
90 tional S given in ([^. In the following section, we find local conservation laws 
corresponding to these variational symmetries. 

3.1. Variational symmetries 

By a variational point symmetry of the action Q , we mean an infinitesimal 
Lie group transformation depending on (t,d,x,p) (no derivatives) 

t ^ t* = t + €^^{t,a,x,p), 
a ^ a* = a + e^{t,a,x,p), 

X ^ X* = X + e fj[t, a, X, p ), 
p^p* =p + eil){t,d,x,p) 

that leaves the action integral (0 unchanged for every region R of spacetime. An 
equivalent requirement for a variational symmetry is for the transformation to 
leave the action differential diS — C At d'®a invariant throughout the integration 
region, such that: 

C{dt>x*,V*x*,p*)At*A^a* = C{dtx,Vx,p)AtA^a. (14) 

A more general definition of variational symmetry allows for C to change by a 
total divergence (see e.g. H): 

C{dt-x*,V*x*,p*)At*A^a* = (£(5tf, Vf,p) -f dtd"a, (15) 
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where = DfN* + DiN^ is a total divergence. 

Extended to act on derivatives of x, this point transformation has infinites¬ 
imal generator X = -\- ipdp + , where 

= Dpjf -x^Dp^'', i=l,...,n, = (16) 

The Lagrangian changes to first order according to £* = £ -f eX£. The volume 
element dt d"a changes to first order in e under the infinitesimal transformation 
according to dt*d"a* = [1-I-dtd^a, which can be shown by direct 
expansion in e or using dSD- Expressing these changes in m, we obtain: 

XC + LD^e = Dpm. (17) 

Einding the variational symmetries of the action integral means solving 
95 invariance condition (II3 for the coefficients r]^, and ip of infinitesimal 

generator X (as well as the unknown fluxes iV*, TV* of the transformation). 

3.2. Solution 

The technique for solving (II3 is to equate all coefficients of derivatives of cc* 
and p to zero, since the point symmetry coefficients i^,rf , and ip only depend 
100 on {t,d,x,p). 

Since the left hand side of dm depends on derivatives of order at most 1 
(i.e. no second or higher derivatives), terms on the right hand side depending 
on higher order derivatives must be independent of the left hand side, so such 
terms must vanish identically in the expression . We conclude that = 

105 N^lt, d,x,p), p = t,l, ... ,n, up to terms that do not influence X. 

We denote C = L + p{J — 1), J = det(a:*). Separating terms in (fT71) that 
are at most linear in x* from those that are at least quadratic in x* yields the 
system 


XL-^ + {L-p)Dp^^^ = Dpm, ( 18 ) 

pXJ + ipj+pJDpC = 0- ( 19 ) 
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From ([B]), XJ = Solving (ITBl) gives 

'^ =-P{il + 'nii), rf = v\t,x), f = f(i,a). (20) 

From the ccja;] term on the left hand side of (IT^ and the ccf and terms on 
the right hand side, 

N* = N; = 0, r=r(a), N^ = N\t,a), i = (21) 

Collecting the remaining terms of zeroth, first, and second order in powers of 
derivatives gives the system 

Piii, - = N* + Ni, 

x\p\=x\Nl, 

xWtTjl, + ^xlxli^j - ^‘) = 0 . 

From the p coefficient in (|2^ and the xlxl coefficients in (|24)) . we get 

^‘(t) = (n + 2)At+ 7 , Cj=nA, vli = \ i = (25) 

where 7 and A are arbitrary constants. The x\xl,i ^ j coefficients in (1241) give 

rf{t,x) = \x'^+uj^^{t)x^+V{t), z = l,...,n, (26) 

where = —w-l®. We find from (|23)) that TV* = a;* 6 ), N"^ = 0, and 6 ®(t) = 
v^t + cr®, where i;®, ct® are constants, and each is defined up to terms that do 
not influence X. 

Finally, for the requirement that = n\ in (I25|) . we can substitute ^® = 
no Aa® + x*(®)- We find that x®(a) is an arbitrary divergenceless vector, or x\ = 0- 
In general, this means that y® = for some anti-symmetric matrix yl®-^ (a) = 
-A^{a). 

In summary, we find the variational point symmetry coefficients of X = 


( 22 ) 

(23) 

(24) 
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+ Tfd^i + tpdp and the fluxes admitted by action functional (O to be: 

^\t) = (n + 2)At + 7, 

C{a) = Xa’’+ x\a), i = 

rf{t, x) = Ax* + + cr*, i = 1,..., n, 

(27) 

i^ip) = -2(n+ l)Ap, 

N\x) = v^x\ 

iV* = 0, f = 1,..., n, 

where A, 7 , 0 ;*-^ = — and cr® are arbitrary constants, and x*(o) is an arbi¬ 

trary vector function with zero divergence: y* = 0 . 

Since each constant is arbitrary, taking the coefficients of these constants 
gives independent symmetry generators. The following symmetries are well 
known in the literature: 

= dt, 

, i 1,..., 7z, 

X^i=td^i, i = (28) 

X^=LO^^X^d,., 

Xoo = X*(d)5ai. 

115 These symmetries respectively correspond to translations in time, translations 
in space, Galilean boosts, rotations in space, and volume-preserving relabelings 
of the particles (the subscript 00 corresponds to “infinite-dimensional”). 

On the other hand, the following point symmetry of action ([9]) and Euler- 
Lagrange equations 0-® has yet to appear in the literature for Lagrangian 
coordinates: 


X\ = {n + 2) t dt + oXdai + x'^d^i — 2{n -I- 1) p dp. (29) 

This symmetry corresponds to a scaling symmetry of the action and Euler- 
Lagrange equations under t —>■ A"'+^t,a* —>• Aa®, X® ^ Ax®, and p —>• A 20+i)p, 
120 A ^ 0. The factors of n appear in the t and p variations due to the volume 
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element dtd"a in the action @ having an n number of a’s. It is known to 
exist for the Eulerian formulation (HD-© as one part of a two-parameter scaling 
symmetry group [^Q], since the action of X on u = Xt is — (n-l- l)u, giving Xa = 
(n -I- 2)t dt + d^i — {n+ 1)m* —2{n+ l)p dp in Eulerian variables. Symmetry 

125 operator Xa is similar to one found by Kambe for an Eulerian variational 
formulation of irrotational compressible flow. Note that (1^^ is also a symmetry 
of the action m for the classical variational formulation of incompressible flow. 


4. Symmetries and Conservation Laws 

In this section we find conservation laws of the Euler equations (01)-® cor- 
130 responding to variational symmetries (P5)) - (l^51) . In addition to the classical 
balance laws of kinetic energy, linear momentum, center of mass, angular mo¬ 
mentum, and the statement of vorticity advection, we find a new conservation 
law for the Euler equations in Lagrangian coordinates corresponding to a scaling 
symmetry of the action functional. 


135 Noether’s Theorem 

We consider a first order Lagrangian C{x^,x^) with independent variables 
and dependent variables x^, fi = t,l,... ,n. The Euler-Lagrange equations 
are given by E^C — 0, where is an Euler operator: 


dxf^ 


-D. 


dxu 




(30) 


Let Xq, = X — -I- ... be the canonical symmetry 

generator, where — x^^’’ is the canonical infinitesimal. The Noether 

operator identity 


2 l|, see also [^or [ 23 , connects Xq, to Ep\ 


Na — Ep + DpRa^p^ (31) 

where Ra^p = -|- .... For variational symmetry Xq, applying (I3ip to C 

in 03 gives 


Dp{R^,pC + eC-N>^) = -a^^EpC. (32) 
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On solutions, {E^C = 0 or =), Eq. (15^ associates to any variational symmetry 
Xq a conservation law (Noether Theorem). For £ in dH]), we have a* = t,x^ = p, 
so (15^ takes the form 

Dt (^a^xi + CC - IV*) + A + f £ - A) = 0. (33) 

Of course, on solutions, we have C = L = \xt\^/2, and J = 1. Thus, from (1551) . 
we obtain: 

a^x{ + -x{x{ ^* — iV*^ + Di + -x{x{ C — = 0. (34) 

4-2. Known conservation laws 

For the time translation symmetry Xj in (1551) . we have C* = 1 and rf = ^'^ = 
ijj = = 0. This gives a* = —xl as the canonical infinitesimal. Substituting 

into the conservation law formula (I34|) . we obtain the conservation form of energy 
balance: 

A + A =0. (35) 

This equation can be recast in Eulerian form. Expanding the time derivative 
and using ([8]), we have 

where m was used. For the flux term, note the following operator identity: 

A = xlD^. = D^. ■ xj - (x^Wkxj) = Aj • xj - (Dk(x^^xj)) = D^. ■ xj, 



x^xlt = u^{ul + 


M'Uli) = A ( 2 



where we used m to conclude that x^^ Dk = Dk ■ x^E This gives: 


D, 


(px^^xi'j = (x’ypx'^^xl^ = D^k (px^) . 


Therefore, in vector notation, we obtain the usual kinetic energy balance law: 


Df 


1 , 


+ V 3 ; 


P+\\u\" ) u 


= 0 . 


(36) 


For the spatial translation symmetries , we have 77 ® ^ and 

A = V® = = 0 - For each symmetry, substitution into conservation law 

formula (|34|) gives: 


D, 


(x^) + A (pi*-*) =0, j 


= 1 ,- n. 


(37) 
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Expanding the derivatives and using xf = 0 recovers momentum equation ([S]) 
in differential form. 

For the Galilean boosts , we have f]^’^ = , and 

= TV* = 0. Substitution into (IMl) gives: 

Dt ^tx{ — x-^'j + Di {tpx'‘^^ =0, j = 1,..., n, (38) 

which, upon integration over d"a, determines the motion of the j —th component 
of the flow’s center of mass. In Eulerian form: 

Dt (tu^) + D^i (tu^u^ + tpS^^ — x^u^) =0, j = 1,... ,n. (39) 

For rotational symmetry X^, we have r]^ = = oj^^x^, and = ijj = = 

0, so dal gives: 

Dt + A {pi^^w^^x’^) = 0. (40) 

In Eulerian form: 

Dt (A*xA*) + D^i (w^^x^u^u* +puj^^x^) = 0. (41) 

If n = 2, dividing by arbitrary constant and letting £ = x^u^ — x^u^ 

gives 

Dti + Ai — px^) + A^ {(-v? + px^) = 0. (42) 

If n = 3, we have w*-! = For each arbitrary constant w*, letting = 

be the angular momentum components, we have 

Dd^ + A- - pe^^^x^) =0, k = \,2, 3, (43) 

which expresses the angular momentum conservation laws. 

140 Symmetry gives a relabeling transformation which maps the fluid labels 
a* to a** = a* + eD(a) in such a way that preserves the volume element d"a: 

=0. Newcomb 23| used this symmetry transformation in conjunction with 
Noether’s First theorem, in spite of the presence of an arbitrary function in its 
symmetry generator. He then derived a conserved quantity using a fixed form 
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145 of the arbitrary divergenceless functions in n = 3 dimensions, this gave a 

representation formula for the fluid vorticity vector Q = x u, first obtained 
by Cauchy till- Note that the arbitrary functions ^*(a) depend on not all in¬ 


dependent variables; see m- For arbitrary functions of not all independent 


variables, Rosenhaus 


24| showed that Noether’s theorem does not yield an infi- 


150 nite set of conserved quantities when such arbitrary functions are involved: only 
one conserved quantity per arbitrary function of spatial variables, and a finite 
number of conserved quantities in the case when the arbitrary functions depend 
on time; see also 2^. Note that, according to Noether’s Second Theorem, an 


infinite symmetry with arbitrary functions of all independent variables would 
155 lead, instead of conserved quantities, to the conclusion that the equations of 
motion are underdetermined; see 


241- 


Using a similar approach to that in 


24| . we show that, instead of an in¬ 


finite set of conserved quantities, this symmetry generates a finite number of 


first integrals of the Euler equations (Cauchy’s invariants [l7|), recovering the 
160 representation formula for uj when n = 3. 

In this case, we have = —xj and = ^*- = ip = = 0. Substitution 

into (IMll gives: 

1 






D, 


) c 


= 0 . 


(44) 


If we use that ^{3) = AJ (a) for arbitrary antisymmetric matrix = —A^ 
and note that 

x'Px^A^ = D.ix’lx'pA^) - D,{x'Px>P)A= = D.ix^x^A^) - x^x^'^ 

(by antisymmetry, x’pjA^^ = 0), we get a density without derivatives of 


A {x%x'yA^) - A 


1 


Dtix^x^)A^ + {p--x^x^ 


A- 


= 0 . 


Let us now integrate over R" and suppose that each A^^{d) is an arbitrary 
function with compact support. By Gauss’ Theorem, the integral of the 
term vanishes, and we obtain 


Dt {x%x’l) d"a = 0. 


(45) 
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By the fundamental lemma of the calculus of variations [18|, the arbitrariness 
of each allows us to conclude that the integrand vanishes for all a. 

For example, if n = 2, we can write A^^{a) = —A'^^ = —A(d), so invoking the 
arbitrariness of A gives: 

A (xitX2 - X2tX^) = 0. (46) 

If n = 3, we can write A^^{a) = e'^^^A^{d), so invoking the arbitrariness of each 
A^ gives 

A = 0, fc = l,2,3. (47) 

This verifies that, rather than an infinite set of conserved quantities, (IT41) leads 
to a finite number, which in this case take the form of first integrals of the 
equations of motion. We evaluate these integrals and show how they correspond 
to the vorticity w. 

If n = 2, we use that = u'- = u'^eXj to rewrite (|46|) as 

{^ 2^1 - ^i^i) = uli {xjxl - x^xl) + ul 2 (xlxf - xlxl) = c(a), 

where c(a) does not depend on t. But in two dimensions, the vorticity is per¬ 
pendicular to the plane: uj = wz, with lo = u^i — We conclude that the left 
hand side is w det(x®) = w by ([ 3 ]), and c = oj. Thus, 

a;(a, t) = w(a, 0), (48) 


165 which shows that the vorticity does not depend on time and is transported with 
the particles by the velocity flow. 

Let us now consider the n = 3 case gZl). We obtain the following first 
integral of the equations of motion ([S]): 

e^^^xJlxT = c\d), (49) 


where c® = 

We now solve (HH)) for the vorticity w. We have xjl = uJ' = u'^gXj, so 
multiplying by x^ gives 


pijk^p q m m 


= C X 


P 

i ’ 
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By a well known formula, det(x®), so using (g]) gives: 

= c^xP. 

But = ojP is the p—th component of the vorticity vector. Since |t=o = 

(JP*, we see that c*(a) = a;*(d, 0) is the initial vorticity. In vector notation: 

w(a, t) = a}(d,0) • Va ai(d, t). (50) 

This representation formula states that the vorticity at any time has been ad- 
vected from the initial vorticity by the Jacobian matrix of Lagrangian map x. 
170 Depending on the flow gradient VaJl, this may lead to enhancement of uj (vortex 
stretching). 


4-3. Scale-Invariance Conservation Law 

Consider now a conservation law corresponding to scaling symmetry (1291) of 
action functional ([H|). Here, 77® = a;®, = a®, = (n + 2)t, ij; = —2(n + l)p, and 

7VP = 0. These infinitesimals produce canonical infinitesimal = x^ — a^xj. — 
{n + 2)txl , so substitution into Noether conservation law (IMl) gives: 


Dt 


x^ - a'=a:i - 




A; 


px'’^ (^x^ — (n + 2)tx{^ 


+ ( 2^14 - P} 


(51) 


= 0 . 


To the author’s knowledge, this conservation law has yet to appear in the lit- 




erature, even in |14| . where the authors performed an exhaustive classification of 


H 


175 the set of conservation laws of the incompressible Euler equations. Kambe 
obtained a similar conservation law for the distinct case of irrotational compress¬ 
ible flow using Noether’s theorem and a scaling symmetry. The main difference 
between his and our conserved densities seems to be the term —a^xl.xl, which 
has no counterpart in his Eulerian formulation. 

Let 


/ = 


7 k 7 

X'^ — a xi — 


■txi 


t > 0. 


(52) 
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be the conserved (constant) integral of conservation law (ISTl) . For conservation 
law (ESI), we let 


E = 


-XfXf 
2 ‘ * 


t > 0. 


(53) 


180 be its conserved integral. Here, E corresponds to the total (kinetic) energy of 
the fluid. 

We show that conserved integral I relates energy E to a, new integral quan¬ 
tity. Equality (l52ll must hold at t = 0. At this time, we have ai(a, 0) = a, so 
, and xl. = 6^^. Therefore, — a’^xj. = 0, so constant I in (l5^ must be 
zero at t = 0, hence for alH > 0. This gives the relation 

J (^x^ - a^xi - t x{^ xl d"a = 0, t > 0. (54) 

Observe that the last term in (1541) can be written as 

where E is as in (j53|) . If we solve ((54ll for E and recall that u = Xt, we find that 

A=Vax) -ud^a, t > 0, (55) 

which represents E as a different integral. 

Although E = / d"a = / d"x is a well defined Eulerian quantity 
(from (IH), d"a = d"x), the integrand in (1551) has no Eulerian counterpart, due 
185 to the term —a ■ VqX • u. Because it is not possible to convert to Eulerian form, 
this integrand is exclusively defined in Lagrangian coordinates. It is interesting 
how an Eulerian quantity can also have a strictly Lagrangian representation. 

We note that (l55)l has a definite limit as t —>■ 0. Let uo(a) = u{a,t) be the 
initial velocity. By (EJ 


1 


^ 1 r 

— 

X — a ■ VaX 

■ u = — 

t 


t 1 


d+tuo-a- Va(a -|- tua) + O(t^) • [mq + 0{t)] 
= (uo -d-VaUo) ■ Uo + 0(t), 


SO taking the limit as t —>■ 0 in (1551) gives 

1 


E = 


n + 2 


uq — a ■ VaUo ) • uq d"a. 


(56) 
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On a related note, since (15^ holds at t = 0, we have E = ^ J \uo\^ d^a. Equating 
this to the right hand side of (|56ll shows that equality holds provided that 


/ 


Va 



d”a = 0. 


(57) 


is satisfied. This is a necessary condition for I to be conserved. In general, it 
will be satisfied when E is finite. 

We now give a physical interpretation for what integral formula (l55l) suggests. 
Since E is conserved, the right hand side of (l55|) must be constant for all time 
t > 0. Therefore, the integral must be proportional to t, which means it becomes 
large as t ^ c». But by Schwarz’s inequality, 


X — a ■ Va x] ■ u d^o 


< 


X — a ■ Va X 


-I 1/2 




so the only part of the integrand in ([5^ that can grow large is x — a ■ 'VxX. 
Specifically, 


X — a ■ Va X 


d"a > 


(n + 2)^ 


Et^ 


t > 0. 


(58) 


190 Since a-Va = is a radial derivative, we interpret x-a-VaX as the radial 


deformation of the particle trajectory x{d,t). The left hand side of (1551) then 
measures the total radial deformation of the fluid at time t. We illustrate this 
with an example. 


Let us consider a well known 
([T|)-(l2]) for n = 2: 


H 


exact solution of the Eulerian formulation 


u^{x, t) = —n(|a;|) x^, 

u^{x,t) = n{\x\) x^, (59) 

p{x,t)=po+ rn(r)^dr, 
do 

where fl(|a;|) is an arbitrary smooth radial function, and po is a constant. Energy 
dSl) takes the form 


pOC 

E = Tr r^n(r)^ dr. 
do 


(60) 
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Solution (15^ describes an ideal fluid with an angular velocity that depends only 
195 upon distance from the origin (a radial eddy). There is no radial velocity, so 
particles transported by the flow retain their original distance to the origin. 

To And the particle trajectories in Lagrangian coordinates, we integrate the 
system xt = u{x, t) and impose the initial condition x{a, 0) = a. Letting 6[a) = 
arctan(a^/a^) be the angle in the Cartesian plane, we obtain 


x^{a,t) = |a| cos (n(|a|)t + 6{a )), 
x^{a, t) = |a| sin (ri(|a|)t + 9{d)). 


(61) 


Observe that 0(|a|) is the angular velocity of the particle given by initial position 
d. For this solution, a direct calculation shows that the radial deformation is as 
follows: 


x^ — d - VaX^ = t |a| 0'(|a|) x^, 
x^ — d ■ VaX? = —t |a| fl^(|a|) x^. 


(62) 


Comparing with (1551) . we see that x — d - VaX = —t|a|n'(|a|) u /0(|a|). 

If O' = 0, or if 0(|d|) is a constant angular velocity, then (l62l) shows that 
the radial deformation is zero. Indeed, in this case, particles of large radius 
200 |d| —>■ oo revolve around the origin as quickly as those near the origin, |d| ^0. 

This represents rigid-body motion, for which no fluid deformation occurs. On 
the other hand, if O' ^ 0, then (1621) shows that radial deformation results. 
Indeed, since the angular velocity O varies with radius |d|, near particles move 
differently from far particles, so the motion is not rigid. In particular, if 0(0) >0 
205 but O —>■ 0 as |d| —>■ oo, then the angular velocity can diminish with distance 
from the origin. In general, this is case for finite energy solutions. 

From (l62l) . there is no fluid deformation for small times t ^ 0, while the 
particles are still at their initial positions. On the other hand, as t increases, the 
fluid deformation increases without bound. Initially close particles of slightly 
210 different radii |a| and |a| -|- <5 become separated from each other due to the 
disparity in their angular velocities 0(|a|) and 0(|a|)-|-5r2'(|a|). The deformation 
increases in (l58l) with energy E because larger E in (l60l) means larger fl, or larger 
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angular velocities. Since 17 ^ 0 as |a| —)• oo, larger fi near |a| = 0 means the 
angles of near and far particles separate more quickly. 

In general, if x-a-VaX = 0, then trajectories are homogeneous in |a|, or a; = 
|a| y{a/\a\,t), where y depends only on direction a/|a| and time t. For example, 
x(a, t) = a (constant flow with zero velocity), or more generally x{d, t) = M{t) a 
for some n x n matrix M(t) with determinant equal to 1 (rigid-body motion). 
Such trajectories may exist for some a in K”, but not all, if Et^ > 0. Moreover, 
since the right hand side of (155)) increases with time, this inequality reflects that 
such trajectories become less important to the overall behavior of the flow as 
time increases. 

Surprisingly, (1581) invalidates the existence of smooth, time-periodic solutions 
of nonzero, finite energy. Indeed, if x{d,T) = x(d,0) = a for some T > 0 and 
all a, then the integral in (1581) vanishes, implying the right side is zero. Since 
r > 0, we conclude that E = 0. Note that m is not a periodic solution if 
it has finite energy (despite being composed of periodic trajectories), since the 
angular velocity changes continuously with a (if 17 = constant, such that x is 
periodic, then E = oo from (l60l) l. Note that there exist time-periodic solutions 
with finite energy for the Navier-Stokes equations (Eulerian formulation) with 


a periodic forcing term 
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27|| . For compressible flow, the existence of time- 


periodic solutions is an open problem 


28|. 


More is true: we cannot have x{a, t + T) = x{a, t) for any t > 0 or T > 0 
(and all a). To show this, we let r > 0,r > 0, and change variables (t,a) —>■ 
(s,6), x{d,t) = y{b,s), such that x(a,T) becomes the initial position of the 
particle trajectories: s = t — T,b = x(a, r), and y{b, 0) = b. But the equations 
of motion (|4|)-(l5]) are invariant under time translations and transformations 
that preserve volume element d"a (cf. X.^ and Xqo in (1551) 1. of which, by (|1]), 
a — ^ x{a,T) is a member. This means y(a,t) is also a solution, so we must 
have y{d, T) ^ y{d, 0) for some a. Since a —>■ x{d, r) is a bijective map (for all 
sufficiently small times; see Lemma 4.4 and Theorem 4.2 in [ll|), this means 
that y{x{di,T), T) ^ y{x(di,T), 0) for some ai. However, by definition of y, this 
means that x{di, t + T) ^ x{di, t), which is what we wanted to show. 
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For this reason, we find that radial deformation of an ideal fluid flow with 
245 finite energy is not time-reversible, such that the bulk flow cannot return to any 
of its previous states. 


5. Conclusion 


We obtained the following results for the incompressible Euler equations in 
Lagrangian coordinates: 

We discussed an action functional dH]) for the incompressible Euler equations. 
Our formulation is a modification of that in [l^ used for compressible flows and 
based off an indirect classical formulation OT. Using this action functional, 
we clarified the relationship between symmetries and conservation laws of the 
incompressible Euler equations; we showed that point symmetries of the action 
result via Noether’s theorem in conservation laws of energy, momentum, center 
of mass, angular momentum, and the statement of vorticity advection. 

It is known H that, for variational systems, there is a one-to-one correspon¬ 
dence between variational symmetries and local conservation laws; therefore, all 
local conservation laws of a given order are generated by corresponding varia¬ 
tional symmetries of this order. In this paper, we have found all lower order 
conservation laws corresponding to point symmetries. However, the conserva¬ 
tion of vorticity, helicity Q], and generalized linear momentum indicate 
the possible existence of generalized (higher) variational symmetries of the Eu¬ 
ler equations in Lagrangian coordinates whose coefficients depend on 

derivatives of x^^p. 

We found a new point symmetry (l29l) of action corresponding to a scaling 
invariance of the action functional and equations of motion. This symmetry is 
similar to one found in [l^ for compressible potential flows. 

Using this symmetry, we constructed a new conservation law m in La¬ 
grangian coordinates. The conserved integral (l52l) relates the total kinetic en¬ 
ergy (l53|l of the fluid, an Eulerian quantity, to a new integral quantity (1551) 
defined exclusively in Lagrangian coordinates. This relationship implies an in- 
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equality (1551) which shows that the total radial deformation of the fluid increases 
with both time and total kinetic energy. The deformation was shown to be time- 
275 irreversible, proving the non-existence of time-periodic solntions with nonzero, 
finite energy. 
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